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Abstract 

New physics effects through the direct CP violation and the decay rate change are 
investigated in the semileptonic D/4 decays, KtiI^u, by including a scalar- 

exchange interaction with a complex coupling. In the decay process, we included 
various excited states as intermediate states decaying to the final hadrons, i^+vr, and 
found that among the intermediate states only the lowest state (K*) is dominant and 
the other higher excited states are negligible, contrary to the B14 decays. We also 
obtained constraints on the new complex coupling within the multi-Higgs doublet 
model and the scalar leptoquark models. 



PACS numbers: 13.20.Fc, 12.60.-i, ll.30.Er 



"'^e-mail: kim@cskini. yonsei. ac.kr, 
^e-mail: jilee@theory. yonsei. ac.kr 
•^e-mail: ngw@cakra. dongguk. ac.kr 



http : / /phya. yonsei . ac . kr / ~ cskim / 



1 Introduction 



Recently we studied the possibility of probing the direct CP violation in the decays, 
Dnl^u |T| and titiI^v 0, in extensions of the Standard Model (SM) by 

including a scalar-exchange interaction with a complex coupling in the weak charged 
current. There we considered, as specific models, the multi-Higgs doublet (MHD) model 
and the scalar-leptoquark (SLQ) models. And we investigated mainly the direct CP 
violation effects, i.e., CP-odd asymmetries for maximally-allowed values of the imaginary 
part of the additional scalar coupling, even though we also found that such a scalar 
coupling could induce a sizable effect to the CP-even total decay rate. In the present 
paper, we take a more general approach and investigate new physics effects by considering 
changes in the (CP-even) total decay rate as well as the (CP-odd) CP violation effects in 
Dn decays, Knl^u. Here we find constraints to both real and imaginary parts of 

the new scalar coupling, by comparing our predictions of the decay rate and the direct 
CP violation effects with the observable experimental results. 

As is well known, in order to observe the direct CP violation effects, there should 
exist interferences not only through weak CP-violating phases but also with different 
CP-conserving strong phases. In the decay Knl^u, we consider it as a two-step 

process: D Ki)lv Knlu, where Ki stands for an intermediate state which decays 

to + vr. In this picture CP-conserving phases come from the absorptive parts of the 
intermediate resonances. The relevant resonance states are K*, Kq and K2 mesons, which 
decay dominantly to Ktt mode (see Table. 1). As shown in Refs. H for the B^ decays, 
the inclusion of higher excited states could amplify the CP violation effects. However, 
we now find that this is not true anymore for the Dn decays, because in the decays 
final state phase spaces are much smaller for those higher excited states than the available 
phase spaces in the -B/4 decays. Therefore, we anticipate that the effects of higher excited 
resonances will be correspondingly reduced and the dominant contributions will come 
mainly through the lowest state K*{892). We will discuss more on this later. 

In Section 2, we present in detail our formalism for the —>■ Knl^u decays within 
the SM and in the models beyond the SM. Section 3 is devoted to the numerical analyses, 
and concluding remarks are also in Section 3. 
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Table 1. Properties and branching ratios of Ktt resonances 



Label i 






nii (MeV) 


Ti (MeV) 









= K*{U30) 


0+ 


1429 


287 


0.62 


1 


Ki 


= K*{892) 


1- 


896 


51 


0.67 


2 


K2 


= k;Iu30) 


2+ 


1432 


109 


0.33 



2 Theoretical Details of the D14 Decays 

We first describe the formalism within the SM for the decay — > Knl^u. Its extensions 
to the models beyond the SM are obtained by rather simple appropriate modifications. 
The decay amplitudes for the processes of Fig. 1 

D~{pn) ^ K^{p^, X^) + W*{q) ^ K+{pk) + 7T-ip^) + l-{pi, Xi) + u{p,) (1) 

are expressed as 

xU,{sJ{K+{Pk)7^-{p^)\\K,{p„X,)), (2) 

where Aj = for spin states (Kq), Aj = ±1,0 for spin 1 states (K*), Xi = ±2, ±1, for 
spin 2 states {K2), and A; is the lepton helicity, ±|. 

The leptonic and hadronic currents are defined, respectively, as 

f = ^,7^(1 -75M, 

= V^e7^(l-75)V^.. (3) 

We assume that the resonance contributions of the intermediate states can be treated by 
the Breit-Wigner form, which is written in the narrow width approximation as 



^^isJ = ^ , ^ , (4) 

where Sj^^ = {px+Pn)'^ and the mj's and Fj's are the masses and widths of the resonances, 
respectively (see Table 1). For the decay parts of the resonances we use 



{K+{PK)7l-{p.)\mp,,x^)) = .jBn,Y,%,,{e*,<p*), (5) 
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D- 




Figure 1: Diagrams for D — > KjW* — > X+tt I vi decays within the SM. 

where Y(^{9, 0) are the J = I spherical harmonics, and the angles 9* and 0* are those of 
the final state 7r~ specified in the Ki rest frame, as defined in Fig. 2c. The couphngs of Ki 
to Xtt are effectively taken into account by the branching fractions, BTZi{Ki X+tt"). 

In order to obtain the full helicity amplitude of the D~ — >■ K'kI~v decay, we first 
consider the amplitude oi D~ Kil'i/i, denoted as M-x/- 

G 

W2 



We express the matrix elements A4x[ ^^^^ following form: 



(6) 



(7) 



where for the decays D KjW* and W* lu, respectively. 



= ew^.{l-{pi,XiMPu)\f^\0), 



(8) 



in terms of the polarization vectors e^y = e(g, A^^) of the virtual W. These e^^'s satisfy 
the relation 



(9) 



(a) D rest frame (b) W* rest frame (c) rest frame 

Figure 2: The decay KiW* — > (A"+7r")(/z/) viewed from the (a) D^, (b) W* and 

(c) Kj rest frames. 



where the summation is over the hehcities A^^ = ±1, 0, s of the virtual W , with the metric 
V± = Vo = -Vs = 1- 

-A, 



We evaluate the leptonic amplitude L^'^ in the rest frame of the virtual W* (see Fig. 
2b). Using the 2-component spinor formalism with explicit polarization vectors [Q, we 
find 

= 2\J^vd±, Lq = —2^J^vdo, = 0, 

= ±2mivdo, Lq = V2miv{d^ — dJ), = — 2m/f , (10) 

where 



mf 1 ± cos 9i 



1 —, d± = 7= — , and dn = sin6'/. (11] 



Here we show only the sign of as a superscript on L. Note that the amplitudes are 
proportional to the lepton mass mi, and the scalar amplitude L~ vanishes due to angular 
momentum conservation. 

For the D ^ Ki transition through the weak charged current 

= V^' - A^", (12) 

the most general forms of matrix elements are, 



for i^o(0+) states 



(i^*fe)|V;|DM)=0, 

{Ko{Pi)\A^\D{pD)) = u+{q'^){pD + Pi)t, + u^{q^){pD - Pi)t,; 
for K*{1~) states: 

{K*ip,,e^)\V^\D{pD)) = ^9{q'')e^.paeriPD+P^nPD-p^r, 
{K*ip,,e,)\A,\D{pn)) = /(g')e* + a+(g^)(et .p^)(p^+p^)^ 



1/1 

^2\f * 



+a^{q ){el ■ Pd){pd ~ Pi) 



for K;{2+) states: 

{K;{p„e2)\V,\D{pD)) = ih{q^)t^,xptrPDa{PD+P.)\pD-Pir, 



{D;{p„ e2)\A^\D{po)) = Kq')tl^yj, + 6+(g')(eL/3P>^)(PD + P 



I) n 



+b-{q'^)i4apPDPD)iPD - Pi] 



ill 



(13) 



where ei and €2 are the polarization vectors of the spin 1 and spin 2 states, respectively. 
Using the above expressions and the polarization vectors f^, we find the non-zero D 
KiW* amphtudes, 



for 2 = 0, Hi 



Si , 



^0 _ „, ^„2\ VQ+Q- 



D ^mJ 



+ U-{q^)Jq^ 



for i = i*'\ i^^ 



^\s,,q' 



f(q^)(ml-s,,-q')+a+iq')Q+Q^ 



(14) 



Vil = f{q')Tg{q')jQ+Q 



V 



for i = 2, i/^^ 



-'o 



rp±l 
-'-±1 



/(g') +a+(g2)(rri 



a-(?')9' 



(15) 



1 




2^6 












1 


Q+Q- r 



k{q^){ml-s„,-q^)+h+{q^)Q+Q^ 



[k{q^)TKq^)^Q+QA, 
k{q') + h+{q^){ml-s,,)+h,{q^)q\ 
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(16) 



where 



A' 



^-7=2 



X=0,s 



A=±l 



where 



The differential partial width can be expressed as 
dT{D- K+TT-l-jyi) 



1 V- I ,AH2 (g'- 



2mD 



where the 4 body phase space (i$4 is 



= dSf^ ■ dif ■ d cos 6* ■ d cos 9i ■ del)* 



(17) 



Combining all the formulae, we can write the SM helicity amplitudes oi D —>■ K^t: I v 
decays as 



(18) 



(19) 



(20) 



(21) 



So far we have established the SM formahsm for the D/4 decays. Now we extend the 
virtual 1^-exchange part in Fig. 1 by including an additional scalar interaction with the 
complex coupling. Then, the decay amplitudes for D~ — > K'^'k~1~vi can be expressed as 

-rai-{vu\)y{vMmHviM\Js\D-{vD)) 

xni(s^)(ir+(p^)7r-(pOII^.(p.,A.)), 



(22) 



where the scalar currents are 
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(23) 



Here the parameter (, which parameterizes contributions from physics beyond the SM, 
is in general a complex number. By using the Dirac equation for the leptonic current, 
Qnj'^ = niijs, the amplitude can be written as 



)) 



yin,{s,,){K^{pK)'n-{pMK^{p^Ai)). (24) 
where the effective hadronic current Vt^ is defined as 

fi^ = J^ + C^J,. (25) 
mi 

In this case the amplitudes M.^. of D ^ KiW have the same form as the previous SM 
case (0) except for the modification in the hadronic current part due to the additional 
scalar current: 

where stands for the hadronic amplitudes modified by the scalar current Jg- Using 
the equations of motion for c and s quarks, we get within the on-shell approximation 



rric — rris rric + 



(27) 



Later we use the approximation, (p^ — P^) ~ {p^ ~p'r ^ = which has been generally 
assumed in quark model calculations of the form factors. After explicit calculation, we 
find that the additional scalar current modifies only the scalar component of 7i^' : 

= (i + C')^f, 

and = Hi' for A^^ = 0, ±1, (28) 



where 

C = —T^—-,C- (29) 

So far, we constructed a formalism for D^^ decays. Since the initial system is not 
CP self-conjugate, any genuine CP-odd observable can be constructed by considering both 
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the D~ decay and its charge-conjugated decay, and by identifying the CP relations of 
their kinematic distributions. Before constructing possible CP-odd asymmetries explicitly, 
we calculate the decay amplitudes for the charge-conjugated process K^ir^l^ui. 
For the charge-conjugated decays, the amplitudes can be written as 

xIl,{sJ{K~ipK)7c+{p.)\\K,{p,,X,)). (30) 

Similarly to the D~ decay, the leptonic amplitudes Z^' for decay are 

Z+ = L+ = -2^vdQ, Lj = 0, 

= ±2mivdQ, Lq = V2miv{d+ — d_), = —2miv. (31) 

And the transition amplitudes 7i^' for KiW* are given by a simple modification 

of the amplitudes Hx' of the D~ decays: 



nl'^=n'^Jg^-g, h^-h, f^^-f^-X^C}, (32) 

which is expected from the property that vector currents change sign under the charge 
conjugation. 

It is easy to see that if ( is real, the amplitude ( p^ of the decay and ( pOD of the 
decay satisfy the CP transition relation: 

A^{9*, 4>*, 9i) = vcpA'^iO*, -4>\ 9i), (33) 

where 6* and (p* in A''^' are the angles of the final state it~^, while those in A''^' are for vr". 
Then, with a complex phase (, dT/d^^ can be decomposed into a CP-even part S and a 
CP-odd part V: 

^ = -(S + V). (34) 

The CP-even part S and the CP-odd part V can be easily identified by making use of the 
CP relation (|33D between D~ and decay amplitudes, and they are expressed as 

S^^±I1, V^'JLZII, ,35) 

d^A d$4 



where F and F are the decay rates for D~ and , respectively. And we use the same 
kinematic variables {s^^-, q'^,9*,9i} for the dT/d^^ except for the replacement of 0* —4>*, 
as shown in Eq. ([33|). The CP-even S term and the CP-odd V term can be obtained 
from decay probabilities. The CP-even quantity S is 

S = 2C{q\sJJ:, (36) 

with 

^={L^SXr\^K*\' + \{V-)UK.\' + \{T-)U^,\' 
+ 2{L^SX')M^K5^K*{V-y + UK*U*^,{T-y) + 2Re{UK*U},,{V-){T-y) 
+ \Il^s\"\L^SX'-i^ + C)LtS%r 

+ \^K^n\{V^)\' + {LtVXn^ + Cr - 2(L+ W)Re((V+))Re(l + (')] 
+ |ni..ni(T+)r + (LtTXril + Cf- 2(L+TfF2°)Re((T+))Re(l + C')] 

- {LtS%')Re{{V^))Re{C) + {Lj SX'){LtV,^Y,')\l + (f] 

+ 2Im(n;,.n;,.)Im((V+))[(L+S° - LtS'M ' {LtS'X)Re{C)] 

+ 2Re(n^.n^.)[(L+5° - L+5°)yo°Re((T+)) - (L+5°yo°)(L+T°r,°)Re(l + (') 

- (L+5X°)Re((T+))Re(C') + (I^^5%°)(L+Tfr,°)|l + Cf] 

+ 2Im(n^.n;,.)Im((r+))[(L+S° - L+S°)Fo° - (L+^°Fo°)Re(C')] 

+ 2Re(n^.n^*)[Re((V+)(r+)*) - (L+T°r2°)Re((V+)) - (L+rfy2°)Re((V+))Re(C0 

- (L+\/°r°)Re((T+))Re(l + O + (L+K°n°)(^^T°r2°)|l + Cf ] 

- 2Im(nK.n;,.)[Im((l^+)(T+)*) - (L+T°r,°)Im((l-+)) - (L+TfF2°)Im((K+))Re(0 
+ (L+K°n°)M(T+))Re(l + C')], (37) 

and the CP-odd quantity T) is 

P = -2Im(C0^^(g',5M)A, (38) 

with 



A 



im((\/+)){(L+K%°)|nK*r + (L+5°Fo°)Re(n^.n;,. 
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+lm{{T^)){{LtTX')\UK.^\' + {LtSX')M^K5llK*) + (L+K°n°)Re(n^.n;,.)} 
+Re{{V^)){{LtTX')M^K*Il},,) - {LtS'X)lm{UK^U*^,)} 

-Re((r+)){(L+K°n°)im(n^.n;,,) + {Lt SX')M^ksIIk^)} 

+ {LtS'X){LtVX')M^K*Il*j,,) + {LtS'X){LtT!Y,')M^Ks^h^] , (39) 

where 

(V^^)^ Y: L^VXi\ E LtTX2\ (40) 

A=0,±1 A=0,±1 

and the overall function C(g^, Sj^^) is given by 

C(i-,sJ = K,P^^ ^^'-f}^ . (41) 



Note that the CP-odd term is proportional to the imaginary part of the parameter ( and 
the lepton mass. Therefore, in order to investigate CP violation in the decays, we 
have to consider massive leptonic {D^ K-Kfi^u) decays. 

3 Numerical Analyses and Conclusions 

First, we calculate i5^4 decay rates through the most dominant resonance, K*{892), by 
varying values of the scalar coupling and then compare those results with the present 
experimental branching ratio 0, 

Bn[D+ (Z^(892)° ^ K-7v+)i2+u^] = (2.9 ± 0.4)%. (42) 

In Fig. 3 we show the allowed parameter space of the complex scalar coupling. For 
example, using the present experimental result of Eq. (|42|), one can constrain the value 
of the scalar coupling down to 

Id -5.0 

at 2-0" level. In our numerical analyses, we use the so-called ISGW2 form factors in 
D Ki transition amplitudes in Eq. (pTS]). 

As mentioned earlier, study of CP violation effects can give a further constraint to the 
imaginary part of the scalar coupling. We consider the so-called optimal observable. An 
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10.0 
B.O 
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4.0 
2.0 
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-4.0 
-6.0 
-8.0 
-10.0 



a — 



CP vlolctioi 



-10.0-8.0 -6.0 -4.0 -2.0 0.0 2.0 4.0 6.0 8.0 10.0 



Figure 3: The allowed parameter space of the complex scalar coupling ( at each confidence 
level, by comparing the theoretical decay rate with the experimental result, Eq. (^), for 
D+ (K*{892f K--iT+)ij,+iy^ decay. The horizontal shaded region represents the 
attainable 2-a limits on the imaginary part through the optimal CP-odd asymmetry (see 
Table 2). 

appropriate real weight function w(s^^, g^; 6**, 6*;, 0*) is usually employed to separate the 
CP-odd V contribution and enhances its analysis power through the CP-odd quantity, 

{wV) = J [wV] . (43) 

And the analysis power is determined by the parameter, 

e^^^. (44) 

For the analysis power e, the number Nd of the D-mesons needed to observe CP violation 
at 1-0" level is 

1 



From the above relation, we can also deduce the bound on the CP-odd parameter for given 
iV^) at an arbitrary confidence level. Certainly, it is desirable to find the optimal weight 
function with the largest analysis power. It is known ^ that when the CP-odd contri- 
bution to the total rate is relatively small, the optimal weight function is approximately 
given as 



V 

s 



-opt 



(f) 

\ iS)' 



(46) 
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Table 2. Attainable 2-cr limits on the imaginary part of the scalar coupling, (, through 
the optimal observable, with the number of D-mesons, = 10^ 



intermediate states 


Attainable 2- 


-a limits 


K*, K* and K* 


|Im(C)l = 


0.121 


K* only 


IMC)I = 


0.126 



We adopt this optimal weight function in our numerical analyses. 

In Table 2, we show the attainable 2-cr limits on the imaginary part of the scalar 
coupling ( through the optimal observable. In order to estimate effects from higher 
excited resonances, we separately present the result obtained by including only the lowest 
state K* as an intermediate state. We can easily see that in D14 decays the effect of higher 
excited resonances is very small. We note the authors of Ref. ||^ analyzed the possibility 
of probing CP-violation by extracting T-odd angular correlations in D ^ K*[-^ Ktt)Ii' 
decay and found that the effects can be detected in some cases. 

In order to get explicit meaning for our analyses, we now consider specific models 
beyond the SM. As specific extensions of the SM, we consider four types of scalar-exchange 
models which preserve the symmetries of the SM Q; one of them is the multi-Higgs- 
doublet (MHD) model [|l^ and the other three are the scalar-leptoquark (SLQ) models 
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(i) Assuming that all but the lightest of the charged scalars effectively decouple 
from fermions, the effective Lagrangian of the MHD model contributing to the 
decay D KttWi is given at energies considerably low compared to Mh by 

C,,„^ = 2V2GFV,s^\m,XZ*icLSR) + m,YZ*{cnSL)]ilRiyL), (47) 

where X, Y and Z are complex coupling constants which can be expressed in terms 
of the charged Higgs mixing matrix elements. From the effective Lagrangian for the 
MHD model, we obtain scalar coupling Ca/hd 



The present bound on Cmhd 

IC^^I < 0.0029 for /i family. (49) 
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Therefore, the MHD bound is aheady too stringent to use this decay mode for 
further constraining the couphng constant. 



(ii) The effective Lagrangians for the three SLQ models |Tl| contributing to the 
decay D Kirlu are written, after a few Fierz rearrangements, in the form 



SLQ 



c 



SLQ 



SLQ 



2 

* 
2 

2M2 



{slCr){vilIr) + ^{sL(y^''CR){viL(T^ylR) 

[slcr){1'ru^l) + \{^L(y^''cR){rRcy^AL) 



+ /i.e.. 



:sL7MCL)(ri7''z/;<i^) + /i.c. 



(50) 



<l>-i 



where j = 1,2 for I = e,n, respectively, and the coupling constants x^j , ylj and 
Zij are in general complex numbers so that the CP symmetry is violated in the 
scalar-fermion Yukawa interaction terms. Then we find that the derived SLQ model 
couplings are 



SLQ 
II 

SLQ 

III 
SLQ 



4V2GFV^tMl 

y2jy2j 



0. 



(51) 



Although there are at present no direct constraints on the SLQ model CP-odd 



parameters in (0), rough constraints on the parameters can be obtained by the 
assumption |jl3| that \x2j\ ~ \x2j\ and \y2j\ ~ \y2j\- That is to say, the leptoquark 
couplings to quarks and leptons belonging to the same generation are of a similar 



size. Then the experimental upper bounds yield [13 



^SLQ 



\QJ < 6.23, \C,J < 15.56 for /i family. 



(52) 



Therefore using these decays, one could extract much more stringent constraints 
on Cg^Q) cis shown in Table 2. 
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To summarize, we have investigated new physics effects through the semileptonic 
decays, — > Kirn^Un, where we extended the weak charged current by including a 
scalar-exchange interaction with a complex coupling. We found that by comparing the 
theoretical decay rate with the observed experimental value, one can constrain both real 
and imaginary parts of the complex scalar coupling. We also investigated the direct CP 
violation effects, and found that one can constrain further the imaginary part through 
the CP-odd optimal asymmetry. We considered as specific models the multi-Higgs dou- 
blet model and the scalar leptoquark models, and found that one can extract much more 
stringent constraints on the scalar- leptoquark couplings, C^^q) through the decay mode 
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